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When a liquid is supercooled towards the glass transition, its dynamics drastically slows 
down, whereas its static structure remains relatively unchanged. Finding a structural 
signature of the dynamic slowing down is a major challenge, yet it is often too subtle to be 
uncovered. Here we discover the structural signatures for both translational and rotational 
dynamics in monolayers of colloidal ellipsoids by video microscopy experiments and 
computer simulations. The correlation lengths of the dynamic slowest-moving clusters, the 
static glassy clusters, the static local structural entropy and the dynamic heterogeneity follow 
the same power-law divergence, suggesting that the kinetic slowing down is caused by a 
decrease in the structural entropy and an increase in the size of the glassy cluster. Ellipsoids 
with different aspect ratios exhibit single- or double-step glass transitions with distinct 
dynamic heterogeneities. These findings demonstrate that the particle shape anisotropy has 
important effects on the structure and dynamics of the glass. 
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The nature of the glass transition is one of the major 
unsolved mysteries in physics On approaching the glass 
transition, the dynamics becomes not only progressively 
slower but also more spatially heterogeneous^'^. Whether such 
dynamic heterogeneity (DH) is associated with any structural 
change is a central problem and an important theoretical 
assumption in glass transition studies^'"*"^. Two types of 
structure signatures, amorphous order^'"*'^'^ and crystalline 
order^, have been proposed as responsible for the dynamics 
arrest. As an archetypal model of amorphous order, icosahedral 
order was often considered to be crucial for vitrification in 
systems of spheres^" and its percolation was associated with 
the glass transition^^. On the other hand, the ground-state 
crystalline structure, even if avoided, is suggested to be important 
for the glass transition^. Recently, crystalline clusters with low 
structural entropy were shown to be responsible for the dynamic 
arrest in systems of spheres^ Despite the intensive search for 
the structural origin of slow dynamics in glass transitions, little is 
known for systems of anisotropic particles. Here we employ 
colloidal ellipsoids that rarely exhibit local crystalline order or 
icosahedral order studied in systems of spheres. By identif)^ing the 
static glassy particles and determining the local structural 
entropies for anisotropic particles, we discover novel structural 
signatures of DHs in both experiments and simulations, and 
confirmed that the glassy clusters with low structural entropy are 
responsible for dynamic arrest not only in translational degrees of 
freedom^^"^^, but also in rotational degrees of freedom. 

Colloids are considered as outstanding model systems for glass 
transition studies because trajectories of individual particles can 
be directly measured by video microscopy^^"^^. Significant 
insights into the glass transition have been obtained from 
experiments on colloidal spheres^ Anisotropic particles, by 
contrast, are much less frequently used as a model system in such 
studies^ despite the fact that they can better mimic molecules 
with anisotropic shapes or interactions. Monodisperse ellipsoids 
are excellent glass formers in two dimensions (2D) because their 
shape can effectively frustrate crystallization and quasi-long- 
range nematic order 

Here we study the glass transitions in monolayers of colloidal 
ellipsoids. The quasi-2D sample enables us to better image and 
track particle motion. The anisotropic shape of the particles 
makes both translational and rotational motion measurable. We 
used four batches of polymethyl methacrylate (PMMA) ellipsoids 
with aspect ratios P = a/^ = 2.3, 3.5, 6.0 and 9.0, where a and b 
are the semi-major and semi-minor axes, respectively. All the 
ellipsoids were obtained by stretching spheres with a diameter of 
2.0 |im. A monolayer of ellipsoids was dispersed in water between 
two parallel horizontal glass walls. For colloids in 2D, the area 
fraction (j)=nahp plays a similar role as the inverse temperature 
1/T in molecular systems^^'^^. Here p is the number density. 
About 10 area fractions within O.2O<0<O.86 were studied for 
each aspect ratio by video microscopy. During the 3- to 6-h 



measurements at each 0, no drift flow or density change was 
observed. The Brownian motion of ~ 6,000 ellipsoids in the field 
of view was typically recorded at one frame per second. The 
center-of-mass positions and orientations of individual ellipsoids 
were tracked using our image-processing algorithm^^. 

In this paper, we focus on the development of DH, the 
evolution of static structures and their relationships towards 
the glass transition at different aspect ratios. First, we characterize 
the kinetic phase behaviour at different aspect ratios by 
measuring the mode- coupling critical point 0^ ^^d the ideal 
glass transition point 0o from the fits to the structural relaxation 
times as a function of density^'^^'^^. Next, we measure the DHs by 
directly visualizing the spatial distribution of fast translational 
(FT), fast rotational (FR), slow translational (ST) and slow 
rotational (SR) particles and by measuring the four-point 
dynamic correlations^'^^. Finally, we analyse the static 
structures by identif)^ing glassy particles and measuring local 
structural entropy. From these analyses in both translational and 
rotational degrees of freedom, we find the following structural 
signatures of the dynamic slowing down in monolayers of 
ellipsoids: the static glassy (that is, high number of neighbours for 
translational motion and high nematic order for orientational 
motion) clusters correspond to low structural entropy and are 
highly correlated to clusters with the slowest dynamics, thereby 
suggesting that the kinetic slowing down is caused by a decrease 
in the structural entropy and a increase in the size of the glassy 
clusters. This indicates that it is the growth of glassy clusters with 
low structural entropy that dictates the dynamic arrest rather 
than the local crystalline order or particular amorphous order, 
which rarely exists in systems of ellipsoids. 

Result 

Phase behaviour. The equilibrium and non- equilibrium phase 
behaviour of ellipsoids in 2D can be qualitatively described as 
follows (Fig. 1). When the aspect ratio P= 1, monodisperse disks 
only form crystals or polycrystals even under the fastest quench 
to high density in experiment and simulation. For P^ 1, ellipsoids 
resemble circular disks and form 2D crystals with free rotational 
motion owing to the spacing between them (Fig. lb). This rotator 
phase or plastic crystal in 2D has quasi-long-range translational 
order and short-range nematic orientational order. As P^oo, 
ellipsoids would form a nematic phase (Fig. le) with quasi-long- 
range orientational order and short-range translational order^^. 
Ellipsoids with intermediate P tend to form local splay structures 
at high density (Fig. If), so that the system has neither quasi-long- 
range translational order nor quasi-long-range orientational 
order, that is, a glass. Moreover, glasses with pseudonematic 
domains may form at large P (Fig. Id), and glasses without 
domains may form at small P (Fig. Ic). The molecular mode- 
coupling theory^^ predicts similar glasses for ellipsoids in three 
dimensions, including an orientational glass at P>2.5 where 
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Figure 1 | Equilibrium and non-equilibrium phases of hard ellipses in 2D. As the aspect ratio P^a/b increases, ellipses form (a) a crystal, (b) a rotator 
phase, (c) a conventional glass without pseudonematic domains, (d) an orientational glass with pseudonematic domains and (e) a nematic phase, 
(f) Three closely packed ellipses tend to have slightly different orientations. 
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rotational motion becomes glassy while the translational motion 
remains quasi-ergodic^^ Simulations about glass of ellipsoids 
mainly focused on small aspect ratios in 2D The recent 

simulation about ellipses in 2D indicates that the rotational glass 
transition should be at a lower density than the translational one 
for P>3, while at a higher density for P<2 (ref. 34). The 
orientational glass has been observed experimentally in quasi-2D 
systems of colloidal ellipsoids with P = 6 (ref. 22). Here we found 
that the 2D orientational glass region shrinks as P decreases and 
completely vanishes at P^2.5. 

We quantitatively measure the glass transition densities 
from the translational and rotational relaxations characterized 
by the self-intermediate scattering function F^{q^t) = 
(^^^^e^^'^^j^^^'^}^^^^) /N and the orientational correlation func- 
tion Ln{t) = (Ej^iCos[n(0^-(t)-O^-(O)])/iV, respectively^^ Xj{t) 
and 6j{t) are the position and orientation of ellipsoid j at time t, q 
is the scattering vector, n is a positive integer, N is the total 
number of particles and < > denotes the ensemble average. We 
chose q — qm corresponding to the first peak in the structure 
factor at high density and n = 3 in Fig. 2. Different choices of q 



and n yield the same glass transition points (Fig. 3c,f). At high 0, 
both Fs(0 and Ljj) develop a two-step relaxation, a signature 
of glassy dynamics. According to the mode-coupling theory, 
the long-time a-relaxation follows the Kohlrausch law e~ {tlx)^ 
(ref. 27) and the structural relaxation time x follows 

T(<^)«(<^,-0)-^ (1) 

where 0c is the critical glass transition point, and y is obtained 
from fitting the a- and |3 -relaxations^ . Translational Tt and 
rotational Xe obtained from Fig. 2a,b can both be well fitted with 
equation (1); (Fig. 2c). The fitted 0^ = 0j = 0.845 ± 0.01 for 
P = 2.3 indicates a single-step glass transition (supercooled 
liquids glass), while (j)l<(i)l for P = 3.5, 6.0 and 9.0 (Figs 2c 
and 3c,f) corresponds to a double-step glass transition 
(supercooled liquid ^ orientational glass ^ glass). This agrees 
with the observation that ellipsoids form pseudonematic 
domains only at P>2.5 such that they no longer rotate but still 
glide translationally within the domains, a signature of an 
orientational glass^ '^^ 
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Figure 2 | Single-step and double-step glass transitions at different aspect ratios, (a) Fs(c/m,0 and (b) mO for ellipsoids with aspect ratio P = a/b = 2.3. 
(c) The translational relaxation tinne xj (solid symbols) and rotational relaxation time xe (open symbols) fitted with equation (1). (d) xj (solid 
symbols) and Xe (open symbols) fitted with equation (2). Inset: the fitted fragility index B. (e) The measured (/)c (equation (1)), (equation (2)), (/)f 
(equation (3)) and (equation (4)) for translational (solid symbols) and rotational (open symbols) motion all show a single-step glass transition at 
P<2.5 and a double-step transition at P>2.5. Curves are guides to the eye. 
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Figure 3 | Translational and rotational relaxations, (a-c) P = 3.5; (d-f) P = 6.0. (a) Fs(qn^,t) at q^ = ^.9\im-\ (d) Fs(c?m,0 at q^^2.3[im-\ (b,e) UCO- 
Legends in (b,e) are the same as those in (a,d), respectively. (c,f) The translational relaxation time xj (solid symbols) and rotational relaxation 
time Te (open symbols) fitted with equation (1) (solid and dashed straight lines). 



Figure 2d shows that the relaxation time is also satisfactorily 
fitted with the Vogel-Fulcher-Tammann equation 

T((/))(xexp[B0/(0o-0)], (2) 

where 0o is the ideal glass transition point and B is the fragility 
index^'"*. Both the fitted 0o ^^d 0^ shown in Fig. 2e confirm the 
single-step glass transition at P<2.5 and the double-step glass 
transition at P>2.5. 0o>0c because all structural relaxations 
cease at 0o, while only free- diffusion relaxations become 
kinetically arrested at 0c (refs 1,3). A smaller B represents a 
more fragile liquid whose viscosity or x is more sensitive to a 
change in T or 0. Fragility is the key concept in the search for a 
universal description of dynamic arrest in glass-forming liquids, 
but its physical origin is still lacking. Recent experiments by 
Mattsson et alP have shown that softer colloidal spheres lead to a 
stronger (less fragile) glass. Here we observed that a larger aspect 
ratio similarly gives a less-fragile glass (inset of Fig. 2d). This 
behaviour was observed in recent simulations on rods^^, but has 
not been tested in experiment before. 

We performed kinetic Monte Carlo (kMC) simulations to 
mimic the Brownian motion of 5,000 ellipses. The same structural 

4 



relaxations, glass transitions and fragility behaviours were 
reproduced as shown in Fig. 4. Figure 4e shows the crossover 
from the single-step glass transition to the double-step glass 
transition at P^l.7, which is lower than the experimental value 
2.5. We attribute this to the following two reasons: (1) an ellipsoid 
carries a layer of water owing to the stick boundary condition, 
making its effective size larger and effective P smaller in 
experiments; (2) surface electric charges on spheres tend to 
redistribute towards the two ends of the resulting ellipsoids, 
which reduces the experimental P. 

DHs. DHs arise from the emergence of domains that are rear- 
ranged in a correlated manner and at different rates It is 
characterized by the cooperatively rearranging regions (CRRs), 
which are closely related to the macroscopic properties of a 
glass CRRs are usually represented by particle clusters with 
FT motion, which have been well studied in systems of 
spheres But much less is known about CRRs that are 

represented by particle clusters with FR^^, ST^^'^^'^^ and SR^"* 
motion. Here we directly visualized the spatial distributions of the 
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Figure 4 | The simulation counterpart of Fig. 2. (a) Fsiq^O and (b) UCt) for ellipsoids with aspect ratio P = a/b = 1.5. (c) The translational relaxation time 
It (solid synnbols) and rotational relaxation tinne xe (open symbols) fitted with equation (1). (d) xjand xe fitted with equation (2). Symbols are the 
same as in (c). Inset: the fitted fragility index B. (e) The measured 0^ 0O/ 0Fand 0s for translational (solid symbols) and rotational (open symbols) motion 
in simulations all show a single-step glass transition at P<1.7 and a double-step transition at P>2.5. 



FT, FR, ST and SR motions of anisotropic particles (Fig. 5a,b,d,e). 
The fastest-moving particles are defined as those with ~ 8% of the 
largest displacements during the cage-breaking time X2 measured 
from the peak of the non- Gaussian parameter of displacements 
Ax: a2(0= <Ax^>/(3<Ax2>2)i^'22'^^. The measured long-time 
diffusion coefficient D((/)) oct2((/>) ~ ^ holds for all densities and 
P's. Hence X2 corresponds to the diffusion time scale^^. The ~ 8% 
corresponds to the non-Gaussian long tails of the displacement 
distribution. Different percentages of fast and slow particles 
yield similar results. The slowest- moving particles are defined as 
those with Ar^(TT) <MSD(4)/4 for translation motion and 
A0^(t0) <MSD(t^)/4 for rotation motion, where t"^ corresponds 
to the middle of the mean-square displacement (MSD) plateau. 
^/MDS{t*) can be regarded as a measurement of cage size^^'^^. 
Fast and slow particles formed the largest clusters at the cage- 
breaking time X2 (ref. 18) and the structural relaxation time x^^, 
respectively. Consequently the fastest-moving particles are 
defined as those with the largest displacements during T2 since 
they dominate diffusion^^, while the slowest-moving particles are 



defined as those with the smallest displacements during x since 
they dominate the structural relaxation^^. For the two samples 
with the highest (/)'s, the slow particles are defined at the longest 
experimental time because their x are longer than the 
experimental time. This leads to be only slightly 

underestimated at the two highest 0's because the longest 
experimental time is on the same order of magnitude as x. 
Subsequently, we can identify the FT, FR, ST and SR clusters by 
using the criterion that two fast (or slow) particles belong to the 
same cluster if they are neighbours. Here we define two ellipsoids 
as neighbours when they overlap after expanding by 1.5 times and 
no other particles intersect at their closest distance^^. 

The FT, FR, ST and SR clusters exhibit different novel 
correlations at P>2.5 and P<2.5. In Fig. 5a,b, we labelled the 
fastest translational particles in green and the fastest rotational 
particles in red. The particles fast in both translational and 
rotational motion are labelled in yellow. In Fig. 5d,e, ellipsoids are 
similarly labelled for slowest particles. The FT, FR, ST and SR 
particles clearly show different positive or negative correlations in 
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Slowest particles 
translation and rotation 
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Figure 5 | Fast, slow and glassy clusters. Spatial distributions of translational (green ~8%) and rotational (red ~8%) fastest (a,b) and slowest 
(d,e) particles, and of translational (g) and rotational (h) glassy (blue ~30%) and slowest (pink ~10%) particles. Particles that are fast (or slow) in both 
translational and rotational motion (a,b,d,e) or that are both slow and glassy (g,h) are labelled in yellow. Scale bar, 20 |im. (a,d), (b,e) and (g,h) are taken 
from the same frame at 0 = 0.81, P = 2.3; 0 = 0.70, P = 6 and 0 = 0.74, P = 6, respectively, (c) The characteristic size of fast clusters (symbols) fitted 
with equation (3) (curves). Inset: the data points collapse onto a master curve with unit slope in the logarithmic plot, (f) The characteristic size 
of slow clusters (symbols) fitted with equation (4) (curves). Inset: the data points collapse onto a master curve with unit slope in the logarithmic plot, 
(i) The characteristic size of glassy clusters (main panel) and the four-point dynamic correlation length ^4 (inset). Both sets of data are well 
fitted with equation (4), except for at P = 2.3. The vertical lines in (c,f,i) mark the onsets of dynamic heterogeneities, below which the data points deviate 
from the fitting curves. 



Fig. 5a,d (P = 2.3) and in Fig. 5b,e (P = 6.0). For P>2.5, the four 
types of clusters are all anticorrelated in space and the 
anticorrelations become stronger as P increases due to better 
developed nematic domains (Fig. 5b,e). For P<2.5, by contrast, 
positive correlations are observed between FT and FR clusters and 
between ST and SR clusters as can be appreciated from the fact 
that most of the clusters are now coloured yellow in Fig. 5a,d. 
Figure 6 is the simulation counterpart of Fig. 5. It shows the same 
results of the spatial distribution and correlations of the FT, FR, 
ST and SR particles. The slow and fast particles in Fig. 6 cover 
about 10% of all particles. The percentage of overlap of fast (or 
slow) transitional and rotational particles (that is, the ratio of 
yellow to (red + yellow) particles) is about 20% (>10%) at 
P= 1.5 and 5% (< 10%) at P = 6.0, which indicates positive and 
negative correlations, respectively. 

The development of DHs can be characterized by the 
increasing sizes of fast and slow clusters. We found that the 
characteristic size (^p of fast clusters is best fitted with^'^^'^^: 

i^{4>) = y(Np) = • ['l>l{4>, - 'l>)]\ (3) 
which diverges at 0f (Figs 5c and 6c); while the characteristic size 

6 



(^s of slow clusters is best fitted with 

= y(Ns) = ^so- </'/(</'s -</'), (4) 
which diverges at 05 (Figs 5f and 6f). Np^s are the numbers of 
particles in a fast cluster and a slow cluster, respectively. 
Equation (4) is consistent with the scaling argument in ref. 4. 
The experimental fitting exponents v = 0.72, 0.59 and 0.448 for ^ 
and v = 0.75, 0.625 and 0.48 for at P = 2.3, 3.5 and 6.0, 
respectively, and the simulation v = 0.38, 0.345 and 0.29 for 
and v = 0.43, 0.385 and 0.33 for at P=1.5, 2.0 and 6.0, 
respectively. These trends indicate that the size of CRRs grows 
faster in more fragile glassy systems. (^f(0) can be better fitted by 
a slightly different form, ^'^^ • (0p — 0) in a broader range of 
(/), but we fit Figs 5c and 6c with equation (4) since it follows 
(;j^ff_ jeff)-v j^CT where Tf^XIcj) for colloids. 0f'^ and 
0g obtained from the best fits of equations (3 and 4) are almost 
identical to and 0q'^, respectively for all P's (Fig. 2e), in 
accordance with the fact that free diffusions dominated by fast 
particles become kinetically arrested at the mode-coupling critical 
density and structural relaxations dominated by slow particles 
cease at the ideal glass transition density 00- Above (jy^^cj)^, the 
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Figure 6 | The simulation counterpart of Fig. 5. The spatial distributions of translational (green) and rotational (red) fast (a,b) and slow (d,e) particles. 
(g,h) Glassy particles (blue) and slow particles (pink) in translational (g) and rotational (h) nnotion. Particles that are fast (or slow) in both 
translational and rotational motion (a,b,d,e) or that are both slow and glassy (g,h) are labelled in yellow. Ellipsoids are plotted as 0.8 times their real size 
(that is, 0.64 times their real area) for clarity. (a,d), (b,e) and (g,h) are taken from the same frame at 0 = 0.85, P = 1.5; 0 = 0.80, P = 6 and 0 = 0.80, 
P = 6, respectively, (c) The characteristic sizes of fast clusters (symbols) fitted with equation (3) in the main text (curves), (f) The characteristic 
sizes of slow clusters (symbols) fitted with equation (4) in the main text (curves), (i) The characteristic sizes of glassy clusters (main panel) and 
the four-point dynamic correlation length (^4 (inset). Both sets of data are well fitted with equation (4), except for the rotational motion at P = 1.5. 
The data points in (c,f,i) collapse onto a master curve in the logarithmic plot in the insets. The solid and dashed vertical lines in (c,f,i) mark the onsets 
of DHs, below which the data points deviate from the fitting curves. 



fast particles dispersed randomly because they were all caged 
separately^^. Moreover, we observe from the insets of Fig. 5f,i that 
the characteristic size of slow clusters (^s reflects the length scale 
of DHs, which is typically characterized by the dynamic 
correlation length (^4 obtained from the four-point density 
correlation function^'^^. The insets of Figs 5i and 6i show that 
(^4 0C 0/(00 — 0)^^ is identical to equation (4) since our ^s^^q- 
Hence (^s(0) (^4(0) for both translational and rotational motion 
and for all P's. 

The structure- dynamics relationship. First, we probe the static 
structure by identifying the glassy particles and show that the 



glassy clusters are highly correlated to the dynamic slow clusters 
in both translational and rotational degrees of freedom. Glassy 
particles are defined as local static structures with strong cage 
effects. We define a particle with Nn>6 nearest neighbours as 
translational glassy, and a particle with a local orientational order 
Sn = cos{2A9j)/N^ > 0.8 as rotational glassy, where AOj is 
the angle difference between the particle and its jth neighbour. 
In practice, the glassy particles in Figs 5g,h and 6g,h are 
defined from the structures averaged over x for better statistics: 

Nn — [/^_^(^^2 ^nl^O^^I/'^T ^ 6 for trauslatioual glassy particles 
and S„ = [f!_{]^, E;=i cosi2A6j{t'))/Nn{n]/Tedt' > 0.8 for 
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rotational glassy particles. Such an average was used in refs 13,14 
and is reasonable since slow particles are defined over the 
structural relaxation time x. Note that and are still static 
quantities after the time averaging. Without time averaging, the 
positive spatial correlations can also be resolved, but with a larger 
noise. and in Figs 5i and 6i are defined without time 
average. Rotational glassy particles have a high local nematic 
orientational order, but the translational glassy particles do not 
exhibit apparent positional order (see the three typical config- 
urations with some splay structures in Fig. 7). 

ATn and exhibited broader probability distributions, that is, 
more structural disorder and stronger frustration, at larger P. This 
observation confirmed the expectation that the lower dynamic 
fragility (inset of Fig. 2d) arises from the stronger structural 
frustration^'^^. The complex local structures in an ellipsoid 
system (Fig. 7) make it difficult to identif)^ whether it is 
geometrical frustration^ or frustration due to competing 
order^'^^'^"*. For geometrical frustration, the slow dynamics 
corresponds to the local stable structure (for example, 
icosahedral order) that cannot tile the whole space, while the 
slow dynamics corresponds to the global stable structure (for 
example, crystalline order) in competing- order frustration. In 
ellipsoids, however, a local favourable structure is ambiguous. The 
close packing of three ellipsoids tends to form a splay structure as 
shown in Fig. If, but close packing of more ellipsoids will form 
different structures (Fig. 7). Hence the observed slow clusters 
(Figs 5d,e and 6d,e) contain various complex structures without a 
particular crystalline or amorphous order. 

The fraction of glassy particles increases with packing fraction 
(j) up to 30% at our highest 0, while the fraction of slow particles 
remains at around 8%. These particles have strong positive 
correlation in space (Figs 5g,h and 6g,h). When P>2.5, 
translational and rotational glassy clusters are spatially antic- 
orrelated: translational glassy particles are mainly found at 
domain boundaries, while rotational glassy particles are mainly 
inside one or multiple adjacent nematic domains. The number of 
particles in a translational (or rotational) glassy cluster is denoted 
as (or ). The mean glassy- cluster size can be well fitted with 

^ \/ (^g) ^^go • 0/(l0ol^) (P^g^ ^^)' except for the 

rotational glassy clusters at P — 23 (Figs 5i and 6i), since this 
system lacks nematic domains (Figs 5a,d and 6a,d) and the glass 
transition is dictated by translational motion. 

The dynamic slowing down is expected to be originated from 
the decrease in the number of accessible configurations^'"*'^. 
Recently, Tanaka et al}^ measured the local structural entropy 
arising from the two-body correlation, S2, in translational degree 
of freedom in systems of hard spheres. They showed that 
crystalline clusters have a lower structural entropy and their sizes 
have the same power-law divergence as the correlation length of 
DH when approaching the glass transition*^. Here we similarly 




Figure 7 | Typical configurations with splay structures, a, b and c exhibit 
some splay structures around a translational glassy particle (in orange) 
with high number of nearest neighbours. 



explore whether glassy clusters are related to low structural 
entropy in systems of ellipsoids. S2 of an ellipsoid / comprises the 
translational and rotational entropies, sn = + > which can be 
evaluated as"*^'"*^ 

poo 

si = - 71/cbP / [gi(r) \ngi{r) - gi{r) + \].rdr, (5) 
Jo 

4i = - ^k^P gMrdr J ^ g.{e\r)\n\giie\r)]d9, (6) 

where is the Boltzmann constant, p is the number density, 
gi{r) is the radial distribution function of centers of mass relative 
to particle iygi{9\r) is the orientational distribution function of the 
angular difference between the long axis of particle / and its 
neighbouring particle at center-of-mass distance r^^. The 
free-volume entropies from multiple layers of neighbours are 
included in equations (5 and 6) via the weights of^/(r) and gi{9\r). 
Both the experiments (Fig. 8a-d) and the simulations (Fig. 9a-d) 
show that the static low-structural-entropy particles and the 
glassy particles are strongly correlated in both translational 
and rotational degrees of freedom. Note that previous 
studies found no correlations between DHs and the local free 
volume in the nearest layer'*^''*'*. By contrast, $2 contains the free- 
volume effect from several layers of neighbours instead of only 
the nearest layer, and remarkably reveals the correlations 
with DHs. 

We measure the spatial correlation of S2, Cs2(r = |rj — rj|) = 
(s2i(rj)s2j(rj)), and fit it with the 2D Ornstein-Zernike 
correlation function Cs2(r) = r" ^/^exp( — r/(^s2)^^'^^ to obtain 
the correlation length ^^2- ^2 in Figs 8e and 9e are averaged over 
|3-relaxation time xp.xp is fitted from the critical- decay law 
F^{q,t) — — hq{t/T^)~^ in the P-relaxation, where the 
nonergodicity parameter and the amplitude /Zq are 
constant^^. Averaging over Xp was also performed when 
defining static crystalline structures for better statistics in 
ref. 46 and was shown to be short enough to avoid possible 
mix up of the static and the dynamic effects. Figure 8e shows 
that (^s2 has the same scaling as Fig. 5f,i and the inset of Fig. 5i 
for all P's, hence ^2(0) (^g((/)) oc (^s(0) (^4(0), with all the 
proportionality prefactors around 1 for both translational 
and rotational motion. Consequently the Vogel-Fulcher- 
Tammann equation x((/)) ocexp[B0/((/)o — 0)] = exp[B42/^s2o] = 
exp[B(^g/(^go] indicates that the kinetic slowing down is related to 
the decrease in the local structural entropy and the increase in the 
size of the glassy clusters. 

Discussion 

Our experiments and simulations show that the correlation 
lengths of the local structural entropy, glassy clusters, slow 
clusters and the four-point dynamic correlation length are all of 
the same order, 42(0)oc^g(0)oc^s(0)oc^4(0)> with the 
same power-law divergence at 0o in both translational and 
orientational degrees of freedom. This set of static- dynamic 
relations and strong spatial correlations of the corresponding 
clusters in Figs 5, 6, 8 and 9 suggest that the static glassy clusters 
((^g(0)) corresponding to the static low structural entropy (^2(0)) 
should be the cause of the slow dynamics ((^s(0)) and the 
increasing DHs ((^4(0)) rather than a series of coincidences. 
This may suggest that the glass transition has a thermodynamic 
origin rather than a purely dynamic origin. The power-law 
divergence of the static correlation length with exponent — 1 
suggests that the glass transition is probably a 2D Ising-type 
critical phenomenon^'^^'^^'"*^'"*^ rather than a random first-order 
transition scenario^, in which the dynamical correlation 
length grows much faster than the static correlation lengths. 
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Figure 8 | Spatial distributions and correlations of local structural entropy, (a-d) The same frame of ellipsoids with aspect ratio P = 6 at (/) = 0.74. 
(a) The spatial distributions of (a) the number of nearest neighbours H^, (b) the local nematic order Sp, (c) the translational structural entropy $2 
and (d) the orientational structural entropy s\. a and c exhibit similar distributions, as do b and d. (e) The correlation leng th of S2'^ fitted with 
<^s2 = (^s2o • 0/(00 - 0)- Only the rotational motion for ellipsoids with P = 2.3 deviates from the fitting curve as in Fig. 5i. (a-d) are not time averaged, 
while (i) is time averaged over the p-relaxation time for lower noises. 



Moreover, the random first-order transition theory concerns 
positional order and not orientational order^, while the scaling 
of our orientational is more consistent with critical 
behaviour. 

In addition, we observed for the first time the crossover from a 
single- step to a double- step glass transition in a system of 
anisotropic particles. The glass transitions for translational and 
rotational motion occur at the same density when the aspect ratio 
P<2.5, but at different densities when P>2.5 where local nematic 
domains start to develop. The translational and rotational CRRs 
exhibit different correlations at P<2.5 and P>2.5. Moreover, 
ellipsoids with a larger aspect ratio lead to stronger frustration and 
a reduced fragility. All the experimental results are confirmed by 
kMC simulations of hard ellipses in 2D. These results show that 



the shape anisotropy of the particles has important effects on the 
phase behaviour and dynamics of the resulting glasses. 

A monolayer of uniformly sized ellipsoids is arguably the 
simplest glass. Note that uniformly sized hard disks can only form 
crystals in 2D. This simple yet versatile glass-forming system with 
measurable translational and rotational dynamics at the single- 
particle level serves as a valuable test bed on which new insights 
into glass transitions in molecular systems can be derived. The 
observed glassy structure and the corresponding CRRs at 
different aspect ratios cast new light on some novel phenomena 
such as photomechanical effects, shape memory and excess 
scattering in molecular, polymeric and liquid- crystal glassy 
systems"*^"^^, where local orientationally ordered domains play 
an important role. Our results pave the way towards a deeper 
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Figure 9 | The simulation counterpart of Fig. 8. (a-d) The same frame of ellipsoids with aspect ratio P = 6 at (/) = 0.80. The spatial distributions of 
(a) the number of nearest neighbours Nn, (b) the local nematic order Sp, (c) the translational structural entropy s[ and (d) the orientational 
structural entropy a and c exhibit similar distributions, as do b and d. (e) The correlation length of s^^^ fitted with = ^^520 • 0/(0o - 0)- Only the 
rotational motion for ellipsoids with P = 1.5 deviates from the fitting curve as in Fig. 6i. (a-d) are not time averaged, while (i) is time averaged over 
the p-relaxation time for lower noises. 



understanding of the glass transition in molecular systems. We 
expect to find similar results in molecular glasses composed of 
non-spherical particles. 

Methods 

Ellipsoids fabrication. We fabricated ellipsoids by stretching non-cross-linked 
PMMA spheres (microparticles GmbH, Berlin) in air using a previously described 

10 



method^^'^^ that is upgraded from stretching spheres in hot oil^'^. Briefly, PMMA 
spheres (0.5% by weight) were added to an aqueous solution of polyvinyl alcohol 
(PVA) (12% by weight) in a Petri dish. After water evaporation, the solid PVA film 
was taken out of the Petri dish and stretched at 130 °C, which is above the glass 
transition temperatures of PVA (rg = 85°C) and PMMA {T^= 105 °C). The 
PMMA spheres embedded in the PVA film were thus stretched into ellipsoids. 
After cooling to room temperature, the PVA was dissolved and washed away with 
deionized water. The cleaned ellipsoid suspension was stabilized with 7 mM 
sodium dodecyl sulphate. The Debye layer was thin and the surface charges had 
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minor effects. From the 2D Gaussian fitting of the images of > 200 isolated 
ellipsoids in a dilute monolayer, we obtained the semi-major axis a and semi-minor 
axis h as fitting parameters. We made four batches of ellipsoids with aspect ratios 
P = a/b = 1.76 |im/0.77 ^im = 2.3, P = 2.33 |am/0.67 ^im = 3.5, P = 3.33 |im/ 
0.56 |im = 6.0 and P = 4.39 |im/0.49 |im = 9.0. The polydispersity of the aspect 
ratio was 5-6% for all P's. 



Sample cells. A monolayer of ellipsoids was dispersed in a thin sample cell 
composed of two parallel glass plates. Glass surfaces of the sample cell were 
rigorously cleaned in a 1:4 mixture of hydrogen peroxide and sulphuric acid by 
sonication. Then the glass was thoroughly rinsed in deionized water. To prepare a 
2D sample, a 0.5- |il droplet of the suspension was placed on a glass slide and then a 
coverslip was placed over it. The droplet spread into a thin film by capillary forces. 
The ellipsoids did not stick to the surfaces. The sample cell was sealed with epoxy 
glue. The area fraction was set by holding the cells up vertically for several hours to 
1 day so that particles slowly drifted downwards under gravity. When the desired 
area fraction was reached in the ~ mm^ chosen area, the cells were laid flat on the 
microscope stage for 2-3 h for equilibration. Two-hour and 24-h equilibration 
resulted in almost the same area fraction, thus the systems were in metastable 
equilibrium. We monitored the central ~ (0.3 mm)^ area for 3-6 h at each density. 
The area fraction remained constant during the whole observation period. For each 
aspect ratio, we measured 3-5 densities in one sample cell and a total of 12 
different densities in three cells. The wall separation was a constant for different 
densities in one cell because we chose a fixed observation area. The wall separations 
in different cells were almost the same since dilute ellipsoids have the same dif- 
fusion coefficients. Under the strong confinement, the heavy ellipsoids always 
stayed in the focal plane and the fluctuation in the z direction was very weak. The 
tips of ellipsoids occasionally overlapped, which shifts the glass transition points to 
< 2% area fraction higher. For brevity, we use the term 2D instead of quasi-2D in 
the main text. The center-of-mass positions and the orientations of individual 
ellipsoids were tracked using our image-processing algorithm^^. The angular 
resolution was 1°, and the spatial resolutions were 0.12 and 0.04 |im along the long 
and short axes, respectively. We projected each step of displacement measured in 
lab frame to the body frame, then connected all steps into a trajectory along the 
long and short axes^ so that transverse motion and longitudinal motion can be 
separately tracked. 

^4{t). We calculated the dynamic correlation length ^^{t) from the four-point 
density correlation function as follows. The four-point correlation function 
of the overlapping particles is defined as g4{r, t) = ^ {J2ij <5(r — rj(0) + r/(0))- 
w{ I rf(0) — ri{t) I) • w(| rj(0) — rj{t) |)), where p is the number density. The 
overlapping function w\ri — is unity if |ri — <0.2a for translation and 
l^^i — 62] <0.l7r for rotation and zero otherwise. We calculated the four-point 
structure factor of overlapping particles as the Fourier transform of g4{r,t): 
S4(q,f) = Jg4{r,t)exp[ — iq ■ r]dr. To extract ^4 at f = x, we fitted S4 with 
54(^,1) = So/[l + {^4q)^]^^. The fitted ^48 are shown in the insets of Figs 5i and 6i. 

Simulation. We performed kMC simulations to mimic the Brownian motion of 
5,000 ellipses. In each kMC step, we randomly generated three types of trial 
moves: translational displacement along the long axis 3^, translational displacement 
along the short axis (5b and rotational displacement dg. follows the uniform 
random distribution in [ — Aa,Aa] . and Sq were similarly generated. The 
maximum displacements A satisfy 



A : A^ A' 



Dh De 



(7) 



where P{6^, P{^h) and P{do) are the acceptance ratios for the displacements ^a> 
and do, respectively. We employed the algorithm in ref. 55 to check the extent of 
overlap of two hard ellipsoids and obtained the acceptance ratios from the 
simulation. The ratios between the short-time diffusion coefficients Da, Db and Dq 
were adopted from the experiments. We also performed simulations with different 
ratios and further confirmed that the short-time dynamics of individual particles 
does not affect the dynamics of the system close to the glass transition points. After 
IVmc MC steps, we have IVmc • ^a(^a) = 2D^At and D = (Da + Db)/2, thus the 
duration of each MC step Af in the unit of the Brownian time Tb can be derived as 



Af/iB = AtD/h" = Nmc • 1 + 



4^2 ' 



(8) 



where Tb is the time taken for an ellipsoid to diffuse a distance h in an infinitely 
dilute suspension. 
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